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Convergence  zone  (CZ)  prediction  models  are  developed  based  on 
acoustic  ray  tracing  theory  as  applied  to  linearly  segmented  sound 
velocity  profiles  (SVP) . The  models  were  developed  into  three 
calculator  programs,  two  for  CZ  range  predictions  under  different 
source  and  receiver  depth  conditions  and  one  for  CZ  gain  and  trans- 
mission loss  (TL)  predictions.  The  performance  of  the  models  as 
programmed  on  Hewlitt-Packard  HP-67  or  HP-97  programmable  calcula- 
tors was  compared  to  the  Fast  Asymptotic  Coherent  Transmission  IF  ACT 
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model  which  is  based  on  similar  but  more  elaborate  theory  and  which 
is  designed  for  use  on  large  digital  computers.  Agreement  of  the 
calculator  programs  with  the  FACT  model  is  fairly  good  when  condi- 
tions are  within  the  design  limitations  of  the  programs  and  environ 
mental  conditions  are  not  unusual. 
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ABSTRACT 


Convergence  zone  (CZ)  prediction  models  are  developed 
based  on  acoustic  ray  tracing  theory  as  applied  to  linearly 
segmented  sound  velocity  profiles  (SVP) . The  models  were 
developed  into  three  calculator  programs,  two  for  CZ  range 
predictions  under  different  source  and  receiver  depth  con- 
ditions and  one  for  CZ  gain  and  transmission  loss  (TL)  pre- 
dictions. The  performance  of  the  models  as  programmed  on 
Hewlitt-Packard  HP-67  or  HP-97  programmable  calculators  was 
compared  to  the  Fast  Asymptotic  Coherent  Transmission  (FACT) 
model  which  is  based  on  similar  but  more  elaborate  theory 
and  which  is  designed  for  use  on  large  digital  computers. 
Agreement  of  the  calculator  programs  with  the  FACT  model  is 
fairly  good  when  conditions  are  within  the  design  limitations 
of  the  programs  and  environmental  conditions  are  not  unusual. 
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a Attenuation  coefficient  (dB/m) . 

C0'C1'C2'* ’ * *ci  Sound  velocities  at  various  points  on  a 

linearly  segmented  sound  velocity 
profile. 


C_  Sound  velocity  at  the  receiver  depth. 


Sound  velocity  at  the  sound  source  depth. 


CZ 


Convergence  zone. 


CZW 


VDl'D2,*’*'Di 


Width  of  the  convergence  zone  annulus. 

Depths  of  various  points  on  a linearly- 
segmented  sound  velocity  profile. 

Depth  of  the  acoustic  receiver. 


Depth  of  the  sound  source. 


DSC 

Ar0,Ar1# ’ * * ,Ari 
^rS'ArR 


Ad 
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Deep  Sound  Channel. 

Horizontal  distances  traveled  by  a sound 
ray  in  traversing  various  layers  of  a 
linearly  segmented  sound  velocity  profile. 

Horizontal  range  corrections  made  to 
sound  ray  cycle  ranges  to  correct  for 
source  and  receiver  depth  separation 
from  the  sonic  layer  depth. 

The  angular  spread  of  all  sound  rays  de- 
parting a sound  source  which  result  in 
convergence  zone  propagation. 

Frequency  (Hz) . 


G Convergence  zone  gain  (dB) . 

gQ,g.,g2, * * ’ ,g.  Sound  velocity  gradients  in  a linearly 

u 1 1 1 segmented  sound  velocity  profile. 


ML 
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Mixed  Layer.  The  upper,  generally  iso- 
thermal, slightly  positive  sound  velocity 
gradient  layer  of  the  ocean. 

The  radius  of  curvature  of  a sound  ray 
in  an  isogradient  layer. 
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Horizontal  range  from  a sound  source  to 
the  inner  edge  of  a convergence  zone 
annulus . 

Horizontal  range  from  a sound  source  to 
the  outer  edge  of  a convergence  zone 
annulus . 

The  cycle  distance  of  a sound  ray  which 
experiences  convergence  zone  refraction. 
It  is  a horizontal  distance  measured 
from  the  point  of  departure  from  the  SLD 
to  the  point  of  return  to  the  SLD. 

The  cycle  distance  for  a ray  which  de- 
parts the  SLD  at  zero  degrees  depres- 
sion angle  from  the  horizontal. 

The  minimum  cycle  distance  of  all  sound 
rays  experiencing  CZ  refraction. 

The  cycle  distance  of  a sound  ray  which 
equals  rQ  but  which  has  a positive 
depression  angle  from  the  horizontal  at 

the  SLD.  The  rays  with  depression 

angles  between  those  for  r,  and  r 
’ 0 rswp 

produce  the  reswept  region  in  a CZ. 

The  working  cycle  distance  used  in  a 
calculator  program  when  iterating  to 
find  a particular  r. 

Sonic  Layer  Depth.  The  depth  of  maximum 
sound  velocity  above  the  DSC  axis. 

Transmission  Loss  (dB) . 

The  TL  in  the  nth  CZ  annulus. 

The  angle  of  departure  of  a sound  ray 
from  the  SLD  measured  downward  from  the 
horizontal . 

The  angle  of  departure  of  the  CZ  sound 
ray  with  minimum  cycle  distance. 
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The  angle  of  departure  of  the  CZ  sound 
ray  which  completes  the  CZ  reswept 
region. 

The  angle  of  arrival  of  a sound  ray  at 
the  receiver  which  departed  the  SLD  at 
zero  depression  angle. 

The  angle  of  arrival  of  a sound  ray  at 
the  receiver  which  departed  the  SLD  at 


rswp 


The  angle  of  departure  of  a sound  ray 
from  the  sound  source  which  reaches 
the  SLD  at  zero  depression  angle. 


The  angle  of  departure  of  a sound  ray 
from  the  sound  source  which  crosses  the 
SLD  at  £ . 


SLD 

when  searching  for  r 


min 


The  average  angle  of  departure  of  the 
sound  rays  within  A6 . 

The  average  angle  of  arrival  of  CZ  sound 
rays  at  the  receiver  depth. 


I.  THE  NEED  FOR  GOOD  CZ  PREDICTIONS 


A.  THE  PROBLEM  OF  THE  AIRBORNE  ASW  UNIT 

In  most  of  the  Pacific  Ocean  and  much  of  the  Atlantic 
Ocean  and  Mediterranean  Sea,  convergence  zone  (CZ)  conditions 
exist  a majority  of  the  time,  and  they  provide  passive  acous- 
tic sensors  with  an  important  means  of  detecting  sounds  emit- 
ted from  submarine  targets.  In  some  areas  the  CZ  regions 
are  the  most  important  contact  regions  considered  in  planning 
acoustic  searches.  Obtaining  accurate  predictions  of  CZ 
sound  propagation  is  therefore  vital  to  the  success  cf  acous- 
tic sensor  tactical  planners. 

Currently,  the  primary  source  of  acoustic  predictions  for 
U.  S.  Navy  units  is  the  Fleet  Numerical  Weather  Center,  Mon- 
terey, California.  Propagation  loss  profiles  for  four  stan- 
dard frequencies  and  three  source  and  receiver  depth 
combinations  are  normally  provided  in  the  ASW  Range  Predic- 
tion System  (ASRAPS)  to  air  ASW  units  when  requested.  The 
profiles,  showing  transmission  loss  (TL)  versus  distance, 
are  generated  on  a large  digital  computer  which  uses  the 
Fast  Asymptotic  Coherent  Transmission  (FACT)  model.  This 
model  uses  as  inputs  a linearly  segmented  sound  velocity 
profile  (SVP) , source  and  receiver  depths,  and  frequencies 
of  interest.  The  SVP  used  can  be  specified  by  the  user  or 
can  come  from  information  stored  at  FNWC  in  the  form  of 
historical  data.  This  stored  data  is  updated  by  bathythermo- 
graph (BT)  reports  through  a complex  weighting  scheme  as  the 


reports  are  received.  Without  going  into  further  detail,  it 
can  be  stated  that  predictions  produced  are  only  as  good  as 
the  data  and  the  computer  model  used,  and  only  as  timely  as 
communications  allow. 

When  entering  a search  area,  a problem  often  arises  con- 
cerning the  TL  profiles  obtained  from  FNWC.  Upon  taking  a 
BT  measurement,  the  unit  often  finds  the  BT  profile  used  to 
generate  the  acoustic  predictions  does  not  agree  with  the 
actual  BT  conditions  in  the  area.  If  this  situation  occurs, 
the  unit  tends  to  lose  confidence  in  the  accuracy  of  the 
predictions  and  tactical  effectiveness  is  felt  to  be  dimin- 
ished by  lack  of  good  information.  The  objective  of  this 
study  was,  therefore,  to  investigate  what  could  be  done  with 
state-of-the-art  programmable  calculators  to  improve  on  the 
in  situ  convergence  zone  predictions  available  to  air  ASW 
units . 

The  reader  only  interested  in  the  calculator  programs 
developed  may  skip  immediately  to  the  appendix. 

B.  ARE  LARGE  COMPUTERS  NECESSARY? 

In  section  5.6  of  Ref.  1,  Principles  of  Underwater  Sound, 
by  R.J.  Urick,  the  author  discusses  the  relative  merits  of 
two  theoretical  approaches  to  obtaining  wave  equation  solu- 
tions in  order  to  describe  the  distribution  of  sound  energy 
in  space  and  time.  Several  references  are  made  to  the  need 
for  digital  computers  to  produce  sound  propagation  descrip- 
tions with  either  theory.  Since  those  comments  were  made, 
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however,  there  haa  been  a revolution  in  the  capabilities  of 
small  programmable  calculators.  Although  it  is  probably  true 
that  computers  are  required  to  produce  a complete  description 
of  sound  propagation  in  the  ocean  with  one  program,  calcula- 
tors are  capable  of  solving  the  different  modes  of  propaga- 
tion one  at  a time  with  separate  programs  to  obtain  a composite 
description.  Examples  of  simple  but  fairly  adequate  calculator 
programs  for  surface  duct,  bottom  bounce,  reliable  acoustic 
path,  and  deep  sound  channel  propagation  modes  are  contained 
in  Refs.  2 and  3.  These  references  also  contain  simple  models 
for  CZ  propagation,  but  they  are  based  on  a mix  of  ray  theory, 
rule-of-thumb,  and  empirical  data.  It  was  felt  that  a better 
CZ  model  needed  to  be  developed. 

C.  DESIRABLE  CHARACTERISITCS  OF  A CALCULATOR  PROGRAM  FOR 

CZ  PREDICTIONS 

1.  The  program  should  require  a minimum  of  easily  avail- 
able input  data.  The  only  information  not  currently  available 
but  which  would  be  needed  by  an  airborne  ASW  unit  is  an  SVP 
from  the  permanent  thermocline  to  the  ocean  bottom  in  the 
unit's  search  area.  A chart  of  this  data  could  easily  be  in- 
cluded in  the  environmental  package  carried  aboard  the  aircraft. 

2.  The  program  should  be  easy  to  operate  and  not  require 
the  operator  to  have  a great  deal  of  insight  into  the  mathe- 
matical model  or  the  internal  operations  of  the  calculator. 

3.  The  output  should  provide  ranges  to  the  inner  and 
outer  edges  of  all  CZ  annuli  of  interest.  It  should  also  pre- 
sent the  expected  TL  for  all  frequencies  of  interest  in  each 
annulus. 
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. The  run  time  of  the  program  should  be  relatively 
short.  This  characteristic  recognizes  that  time  is  impor- 
tant to  an  on-station  unit. 

5.  The  program  should  be  based  on  generally  accepted 
acoustic  theory.  This  characteristic  is  desirable  because  a 
user  would  probably  have  more  confidence  in  such  a model 
than  one  based  on  empirical  data  and  thus  applicable  only  to 

a specific  ocean  basin.  With  empirical  models,  the  user  often 
wonders  if  the  area  he  intends  to  search  corresponds  to  the 
mean  set  of  conditions  used  to  generate  the  model  or  is  some- 
how different. 

6.  Ideally,  the  program's  performance  should  agree 
closely  with  the  generally  accepted  large  computer  models 
currently  in  use. 


II.  USING  ICAPS  AS  THE  STANDARD  FOR  COMPARISON 


A.  RATIONALE  FOR  ONLY  ONE  LOCATION  PER  OCEAN  BASIN 

As  will  be  demonstrated,  the  CZ  characteristics  of  the 
three  locations  studied  vary  considerably.  The  deep  sound 
channels  which  produce  CZ  phenomena  are  quite  different  as 
are  the  ranges  from  source  to  CZ  annuli.  The  objective  of 
this  study  was  to  produce  a mathematical  model  for  CZ  pre- 
dictions for  use  on  small  programmable  calculators.  It  was 
reasoned  that  if  the  model  would  work  for  the  different  con- 
ditions of  the  three  locations  studied-  it  would  work  for  all 
of  the  variations  to  be  expected  within  any  one  of  the  ocean 
areas. 

B.  DESCRIPTION  OF  ICAPS 

The  Integrated  Carrier  ASW  Prediction  System  (ICAPS)  is 
a passive  and  active  acoustic  prediction  system  developed 
for  installation  aboard  aircraft  carriers  and  other  large 
naval  vessels  which  have  digital  computers.  It  contains 
four  sets  of  historical  environmental  data,  one  each  for  the 
North  Pacific,  North  Atlantic,  and  Indian  Oceans,  and  one 
for  the  Mediterranean  and  Black  Seas.  It  also  contains  sev- 
eral production  programs  for  predicting  naval  sonar  system 
performance.  The  FACT  model  is  used  in  the  passive  sensor 
predictions.  This  is  the  same  model  used  at  FNWC  for  ASRAPS. 
Reference  4 contains  a description  of  the  installation  and 
operation  of  ICAPS  in  the  IBM  360  Computer  Center  at  the 
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Naval  Postgraduate  School.  Reference  5 contains  a descrip- 
tion of  the  mathematics  used  in  the  FACT  model. 

C.  DESCRIPTION  OF  HISTORICAL  ENVIRONMENTAL  DATA  FILES 

Figures  1(a)  through  1C1)  depict  twelve  sound  velocity 
profiles  produced  by  ICAPS  from  its  historical  environmental 
data  files.  Figures  1(a)  through  1(d)  show  SVP  information 
for  the  months  of  February,  May,  August,  and  November  for 
40N  140W  in  the  Pacific  Ocean.  Figures  1(e)  through  1(h) 
are  the  same  information  for  3 IN  69W  in  the  Atlantic  Ocean, 
and  likewise,  Figures  l(i)  through  1(1)  are  for  36N  18E  in 
the  Mediterranean  Sea. 

The  historical  data  files  used  to  produce  these  profiles 
consist  of  temperature  and  salinity  values  for  over  thirty 
depths,  for  four  seasons  of  the  year,  and  for  many  locations 
spaced  at  one  to  five  degree  latitude  and  longitude  inter- 
vals in  each  ocean  area  covered.  When  specific  latitude, 
longitude,  and  date  are  specified,  interpolations  are  per- 
formed to  produce  the  approximate  temperature  and  salinity 
profiles  to  be  expected  at  that  location  and  date.  This  in- 
formation is  then  converted  to  an  SVP . The  output  from  this 
portion  of  the  system  consists  of  seven  columns  of  values, 
one  each  for  depth  in  meters  and  feet,  temperature  in  Cel- 
sius and  Fahrenheit,  salinity,  and  sound  velocity  in  meters 
per  second  and  feet  per  second.  The  depths  associated  with 
these  quantities  begin  at  ten  meter  intervals  near  the  ocean 
surface  and  gradually  increase  through  25,  50,  100,  250,  500, 
and  1,000  meter  intervals  as  depth  increases.  The  last  line 
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of  values  is  for  the  ocean  bottom  depth  which  was  part  of 
the  input  data. 

D.  COMPARISON  OF  DEEP  SOUND  CHANNEL  CHARACTERISTICS 

Figure  2 depicts  the  deep  sound  channel  (DSC)  portion 
of  the  May  SVP  for  the  Pacific,  Atlantic  and  Mediterranean 
coordinates  mentioned  earlier  in  a composite  graph  drawn  to 
scale.  As  can  be  seen  in  that  figure,  DSC  characteristics 
of  the  three  areas  differ  considerably.  The  vertical  ex- 
tent of  the  channels  varies  from  1100  meters  in  the  Mediter- 
ranean to  4200  meters  in  the  Atlantic.  The  change  in  sound 
velocity  between  sonic  layer  depth  (SLD)  and  DSC  axis  (point 
of  minimum  velocity)  varies  from  15  m/sec  in  the  Mediterra- 
nean to  38  m/sec  in  the  Atlantic.  The  depth  of  the  DSC  axis 
varies  from  100  meters  in  the  Mediterranean  to  1300  meters 
in  the  Atlantic.  Pacific  Ocean  values  are  between  the  others 
for  all  of  those  characteristics . Sound  velocity  near  the 
surface  is  much  greater  in  the  Atlantic  and  Mediterranean 
than  in  the  Pacific,  and  sound  velocity  near  the  bottom  of 
the  three  basins  (not  shown  in  the  figure)  is  about  6 m/sec 
greater  in  the  Atlantic  than  in  the  Pacific  and  about  52 
m/sec  greater  in  the  Mediterranean  than  in  the  Pacific  at 
equal  depths.  Also  note  the  subsurface  sound  channel  lo- 
cated about  100  to  500  meters  below  the  surface  in  the  At- 
lantic profile. 
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Figure  2.  Comparison  of  Pacific,  Atlantic,  and 
Mediterranean  Deep  Sound  Channel  Characteristics. 
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E.  COMPARISON  OF  CZ  CHARACTERISTICS  IN  THREE  OCEANS 
1 . Method  Used  in  Obtaining  Data  for  Comparison 

There  were  two  primary  objectives  in  gathering 
twelve  ICAPS  runs  from  each  of  the  three  ocean  areas.  First, 
it  was  desired  to  obtain  sufficient  data  to  determine  which 
CZ  characteristics  are  common  to  all  areas  and  which  charac- 
teristics are  peculiar  to  specific  basins.  Secondly,  it  was 
desired  to  obtain  a standard  of  comparison  for  any  calculator 
program  which  might  be  developed.  To  fulfill  the  first  ob- 
jective, it  was  decided  to  keep  the  input  variables  the  same 
in  all  areas,  varying  them  one  at  a time,  in  order  to  better 
compare  the  differences  observed  in  the  various  runs.  For 
each  of  the  three  locations,  the  inputs  varied  were  season 
of  the  year  and  source  and  receiver  depth  combination.  Re- 
ceiver depths  of  60  and  300  feet  and  source  depths  of  60  and 
400  feet  were  used.  Each  of  the  ICAPS  outputs  consisted  of 
TL  profiles  for  four  frequencies  out  to  a range  of  250  kyds . 

Originally,  it  was  intended  to  collect  twelve  data 
from  each  profile.  These  data  were  to  be  the  range,  width, 

CZ  gain,  and  Transmission  Loss  for  each  of  the  first  three 
convergence  zone  annuli.  As  it  turned  out,  somewhat  less 
data  was  collected  and  tabulated.  There  were  several  rea- 
sons for  this.  First,  the  February  SVP  in  the  Mediterranean 
contained  no  sound  channel  and  therefore  no  convergence  zones 
existed.  Secondly,  all  of  the  third  CZ  data  for  the  Atlan- 
tic was  thrown  out  on  the  grounds  that  it  was  almost  always 
the  same  and  that  it  was  inconsistent  with  information  from 
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the  first  two  CZ  annuli  in  any  particular  profile.  The 
reason  for  this  occurrence  is  not  know.  Finally,  it  was 
impossible  to  obtain  some  of  the  desired  data  because  of 
the  smooth  way  in  which  the  CZ  path  blended  with  other  com- 
petitive propagation  modes.  One  could  not  tell  what  was 
CZ  and  what  was  not  in  those  cases. 

2 . CZ  Range  and  Width  Analysis 

The  transmission  loss  profiles  produced  by  ICAPS 
are  presented  in  two  formats,  a table  of  TL  values  for  each 
kiloyard  of  range  from  the  source  and  a graph  of  the  same 
information.  Because  the  TL  values  are  tabulated  at  kilo- 
yard  intervals,  it  is  impossible  to  be  more  accurate  than 
that  interval  in  determining  where  a CZ  begins  and  ends. 
AI30,  it  was  difficult  to  be  consistent  in  picking  the 
points  representing  the  edges  of  CZ  annuli  because  of  the 
variety  of  graph  shapes,  TL  levels,  and  other  propagation 
mode  interferences.  In  any  event,  an  attempt  was  made  to 
satisfy  one  basic  criterion  in  choosing  leading  and  trailing 
edges  of  the  annuli:  Do  these  ranges  best  represent  the 
apparent  location  of  the  annulus  regardless  of  the  TL  levels 
involved?  Admittedly,  the  ranges  picked  were  often  based 
on  subjective  judgement,  and  it  cannot  be  stated  with  com- 
plete certainty  that  only  CZ  mode  propagation  contributed 
to  the  TL  peaks  judged  to  be  the  CZ  annuli. 

Table  I contains  the  CZ  range  and  width  data  that 
could  be  gleaned  from  the  ICAPS  profiles.  In  the  table, 

RCZi  is  the  range  to  the  inner  edge  of  the  first,  second,  or 
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Table  I.  1CAPS  CZ 
Range  A Width  Data 
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third  CZ  annulus  to  the  nearest  one  half  nautical  mile.  CZW 
is  the  width  of  the  annulus  also  the  nearest  one  half  nauti- 
cal mile.  The  third  column  of  numbers  is  the  ratio  of  CZW 
to  the  range  of  the  outer  edge  of  the  annulus  (RCZo) , ex- 
pressed as  a percentage. 

After  carefully  studying  this  data,  the  following 
conclusions  were  made  concerning  CZ  propagation: 

(1)  The  range  to  the  first  CZ  is  approximately 
14  to  18  nai  at  the  Mediterranean  location,  23  to  27  nm  at 

the  Pacific  location,  and  33  to  35  nm  at  the  Atlantic  location. 

(2)  Range  to  the  CZ  decreases  and  annulus 
width  increases  as  source  and  receiver  get  deeper  in  all 
cases . 

(3)  The  ranges  to  the  second  and  third  annuli 
are  approximately  whole  number  multiples  of  the  ranges  to  the 
inner  and  outer  edges  of  the  first  annulus  in  all  cases. 

(4)  The  range  or  width  of  a CZ  annulus  does 
not  appear  to  have  any  significant  frequency  dependence. 

3 . CZ  Gain  and  Transmission  Loss  Analysis 

Convergence  zone  gain  is  defined  as  the  difference 
between  the  transmission  loss  expected  under  conditions  of 
spherical  propagation  and  the  actual  transmission  loss  ob- 
served. This  definition  is  expressed  in  Eq.  (1)  . 

G = 20  log(r)  + a(.r)  - TL  (1) 

In  this  equation,  G is  the  CZ  gain,  r is  the  range  to  the 
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CZ  annulus,  a is  the  attenuation  coefficient  associated  with 
the  frequency  of  interest,  and  TL  is  the  actual  transmission 
loss  observed  in  the  CZ  annulus  for  that  frequency.  All 
terms  in  Eq.  (1)  are  in  decibels  (dB) . 

In  actual  convergence  sones,  TL  (and  therefore  gain) 
is  by  no  means  a constant  value.  Contributions  of  several 
possible  propagation  paths  at  any  one  point  and  the  time 
varying  nature  of  sound  paths  in  the  ocean  cause  coherence 
effects  to  exist.  These  effects  make  TL  vary  in  both  space 
and  time.  Coherence  effects  are  more  pronounced  at  lower 
frequencies  (longer  wavelengths)  where  the  time  varying  ef- 
fects are  small  compared  to  spatially  distributed  effects. 

In  ICAPS,  the  more  predictable  coherence  conditions  are  in- 
cluded in  the  mathematical  model. 

Since  a single  TL  value  was  desired  for  the  envi- 
sioned calculator  model,  an  attempt  was  made  to  pick  the 
"average"  TL  in  the  ICAPS  CZ  annuli.  As  with  the  range  es- 
timates, this  called  for  subjective  judgement.  Figure  3 
shows  a typical  ICAPS  CZ  presentation  which  has  coherence 
effects  in  evidence.  The  figure  suggests  how  an  "average" 

TL  was  chosen  as  best  representing  that  annulus.  'Two  levels 
were  chosen  (labeled  high  and  low  in  the  figure)  which  bracket 
the  majority  of  the  TL  points  within  the  annulus.  The  ap- 
proximate midpoint  between  those  levels  was  then  picked  as 
"the"  TL  for  that  CZ. 

As  an  extra  point  of  interest,  the  high  and  low  TL 
levels  were  studied.  It  was  noted  that  ICAPS  predicts  TL 


Transmission  Loss 


Figure  3.  Estimating  Transmission  Loss  in  a C2 
Annulus  from  an  ICAFS  TL  Frofile. 
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variations  from  about  t lOdB  around  the  "average"  level.  If 
this  is  truly  representative  of  CZ  coherence  effects,  an  ASW 
unit  armed  only  with  an  estimate  of  the  "average"  TL  in  a 
certain  CZ  annulus  should  expect  to  see  variations  of  about 
that  magnitude  around  the  estimate  in  hand. 

Using  TL  levels  estimated  by  the  procedure  described 
above,  and  employing  Eq.  Cl) , CZ  gain  values  predicted  by 
ICAPS  were  obtained  and  tabulated.  The  values  produced  are 
contained  in  Tables  II  and  III.  Table  II  shows  all  of  the 
data  from  the  Pacific  location.  Table  III  contains  only  300 
Hz  data  from  the  Atlantic  and  Mediterranean  locations.  (50, 
850,  and  1700  Hz  data  were  omitted  from  Table  III  because  it 
became  obvious  during  data  collection  that  G is  not  frequency 
dependent. ) 

Again  after  careful  study,  the  following  conclusions 
were  drawn  concerning  CZ  gain: 

(1)  CZ  gain  values  range  from  about  eight  to 
twenty  dB  in  all  three  areas  observed. 

(2)  CZ  gain  is  the  same  value  for  first, 
second,  and  third  CZ  in  any  given  case. 

(3)  In  general,  CZ  gain  is  independent  of 
frequency.  An  exception  to  this  conclusion  is  that  at  low 
frequency  (below  300Hz)  , especially  when  source  or  receiver 
or  both  are  above  the  SLD  and/or  near  the  surface,  there  is 
apparently  somewhat  less  gain  than  evident  for  higher  fre- 
quencies. This  difference  is  probably  due  to  stronger  dif- 
fraction of  the  longer  wavelengths. 
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(4 ) CZ  gain  seems  to  be  highest  when  source 
and  receiver  are  at  or  near  the  same  depth. 
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Month 
&LD ( ft) 


FEB 

tt? 

i 


Revr/Tqt 


60/60 


50 

300 

850 

1700 


CZ  Gain  (dB) 

1st  CZ  2nd  CZ  3rd  CZ 


7 7 6 

14  13  12 

14  13  13 

15  15  12 


60/400 


50 

300 

850 

1700 


9 

12 

11 

11 


9 9 

12  11 

11  12 

11  12 


300/400 


50  12 

300  12 

850  10 

1700  12 


11 

12 

11 

10 


10 

11 

12 

11 


60/60 


50  16 

300  16 

850  lb 

1700  16 


17 

16 

17 

16 

17 

16 

17 

17 

MAY 

TT 


60/400 


50  10 

300  10 

850  11 

1700  9 


10 

11 

12 

11 


12 

13 

15 

12 


I 


300/400 


50 

300 

850 

1700 


14 

13 

13 

12 


14 

14 

13 

14 

14 

11 

11 

13 

Table  II.  ICAPS  CZ  Gain  Data  Observed  at  the  Pacific  Ocean 
Location.  (Page  1 of  2) 
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Month 

Rcvr/Tqt 

Freq 

Thz) 

CZ 

Gain  (dB) 

SLDTTt) 

(it) 

1st  CZ 

2nd  CZ 

'3rd  CZ 

50 

11 

11 

13 

60/60 

300 

850 

14 

15 

13 

18 

12 

20 

1700 

16 

17 

18 

50 

7 

7 

8 

AUG 

0 

60/400 

300 

850 

9 

12 

12 

14 

13 

15 

1700 

12 

13 

15 

50 

14 

14 

14 

300/400 

300 

11 

14 

12 

850 

11 

12 

15 

1700 

12 

13 

11 

50 

8 

7 

6 

60/60 

300 

9 

12 

12 

850 

15 

18 

16 

1700 

14 

15 

15 

50 

9 

9 

8 

NOV 

~TS 

60/400 

300 

850 

12 

11 

13 

11 

11 

11 

1700 

11 

11 

11 

50 

11 

11 

11 

300/400 

300 

850 

12 

9 

11 

11 

11 

12 

1700 

10 

9 

11 

Table  II.  ICAPS  CZ  Gain  Data  Observed  at  the  Pacific  Ocean 
Location.  (Page  2 of  2) 


MEDITERRANEAN  ATLANTIC 


Ocean 


Month 
SLD ( ft) 


FEB 


300  Hz  CZ  Gain 


(ftp 

1st  CZ 

2nd  CZ 

“3rd  CZ 

60/60 

17 

14 

60/400 

13 

13 

300/400 

17 

18 

MAY 

IT* 


60/60 

19 

14 

60/400 

16 

12 

300/400 

13 

17 

AUG 

0 

60/60 

60/400 

300/400 

15 

10 

14 

17 

13 

15 

60/90 

13 

11 

NOV 

60/400 

12 

10 

0 

300/400 

13 

12 

MAY 

60/60 

14 

13 

16 

Iff" 

60/400 

12 

10 

10 

300/400 

17 

18 

20 

AUG 

TT 

60/60 

60/400 

300/400 

17 

11 

19 

9 

17 

10 

NOV 

60/60 

11 

13 

15 

IP 

60/400 

11 

12 

300/400 

12 

14 

15 

Table  III. 


ICAPS  300Hz  CZ  Gain  Data  Observed  at  Atlantic 
and  Mediterranean  Locations. 
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III.  CZ  RAY  THEORY  ANALYSIS  AND  MODEL  DEVELOPMENT 


A.  CZ  RANGE  AND  WIDTH 

It  was  decided  to  use  ray  tracing  as  the  method  for  de- 
termining CZ  range  and  width  because  of  the  simplicity  of 
the  mathematics  involved  and  because  of  the  intuitive  appeal 
of  sound  rays  depicting  the  propagation  of  sound.  The  alter- 
native approach,  that  of  normal  mode  theory,  was  rejected  on 
the  grounds  that  it  would  be  much  more  complicated,  requiring 
capabilities  far  beyond  those  available  in  the  calculators 
at  hand. 

Figure  4 shows  four  sound  rays  of  particular  interest  in 
CZ  propagation.  The  order  of  these  rays  is  described  for  the 
"typical"  case  in  the  following  discussion.  An  "atypical" 
case  will  be  mentioned  later. 

Ray  #1  departs  the  SLD  at  zero  degree  depression  angle. 

It  reaches  its  greatest  depth  at  the  bottom  of  the  DSC  and 
returns  to  the  SLD  at  some  particular  range  and  at  horizontal 
incidence.  The  horizontal  range  from  SLD  to  SLD  is  termed 
cycle  distance.  The  cycle  distance  for  this  ray  is  desig- 
nated rQ. 

Ray  #2  is  the  next  ray  of  interest  found  as  the  departure 
angle  from  the  SLD  is  increased  downward.  This  ray  passes 
down  through  and  below  the  bottom  of  the  DSC  before  turning 
back  upward.  It  returns  to  the  SLD  at  the  shortest  range 
from  the  starting  point  of  any  ray  within  the  bundle  of  rays 


undergoing  CZ  refraction.  This  range  is  designated  r . . 

mm 

The  angle  of  departure  for  this  ray  is  designated  9 . . 

man 

Each  ray  between  #1  and  #2  crosses  all  of  the  previous  (les- 
ser departure  angle)  rays  on  its  way  up  from  its  lowest 
depth. 

As  departure  angle  from  the  SLD  is  further  increased, 
the  next  sound  ray  of  interest,  #3,  is  located.  This  ray 
has  a cycle  distance  equal  to  rQ.  Its  maximum  depth  is 
greater  than  that  for  ray  #2.  It  departs  from  and  arrives 
back  at  the  SLD  at  an  angle  designated  9rsWp*  Rays  between 
#2  and  #3  do  not  cross  each  other,  but  they  do  cross  the 
earlier  rays  on  their  way  back  up  to  the  SLD. 

In  the  CZ  annulus,  the  rays  between  #1  and  #2  sweep  in- 
ward toward  the  source  as  departure  angle  increases.  After 
ray  #2  they  sweep  out  away  from  the  source  as  angle  increases 
further.  For  this  reason,  the  region  formed  by  rays  between 
#1  and  #3  is  called  the  reswept  zone. 

Finally,  as  angle  of  departure  from  the  SLD  is  increased 
to  maximum  angle  for  CZ  propagation,  we  observe  ray  #4.  This 
ray  turns  upward  at  a depth  equal  to  the  water  column  depth 
at  that  location.  It  returns  to  the  SLD  at  the  greatest  dis- 
tance of  all  CZ  refracted  rays.  Rays  departing  the  SLD  at 
angles  greater  than  that  for  ray  #4  would  be  reflected  off 
the  bottom  and  are  not  of  significance  for  the  CZ  propagation 
path. 

As  mentioned  earlier,  this  progression  of  rays  exists 
in  a "typical"  CZ  situation.  If,  however,  the  ocean  bottom 
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were  more  shallow,  cycle  distance  for  ray  #4  would  be  re- 
duced. If  the  bottom  were  shallow  enough,  ray  #4's  cycle 
distance  would  be  less  than  r^.  In  that  case,  the  reswept 
zone  would  be  reduced  to  the  region  between  rays  #2  and  #4. 
This  situation  is  called  the  "atypical"  case. 

B.  RANGE  AND  WIDTH  MODEL 

In  the  calculator  programs  developed,  provision  is  made 
for  entering  and  storing  a five  point  sound  velocity  profile 
which  defines  the  DSC  only.  The  first  depth  and  velocity 
pair  entered  (D^C^)  equate  to  the  appropriate  values  found 
at  the  SLD . The  fifth  depth  entered  (Dg)  is  the  depth  at 
the  bottom  of  the  DSC  where  sound  velocity  is  equal  to  that 
at  the  SLD.  The  other  three  depth/velocity  pairs  must  be 
picked  subjectively  from  a graph  of  the  SVP  of  interest.  If 
a mixed  layer  exists,  the  gradient  in  that  layer  is  taken  to 
be  0.02  sec  ^ (a  purely  pressure  induced  gradient  to  two 
place  accuracy) . The  program  calculates  the  four  layer  gra- 
dients within  the  DSC  profile  entered,  and  uses  the  fourth 
(deepest)  layer  gradient  in  ray  calculations  that  occur  below 
the  DSC.  It  would  have  been  desirable  to  allow  several  more 
points  in  the  SVP,  but  calculator  data  storage  capacity  and 
program  step  limitations  preclude  more  than  five  depth/veloc- 
ity pairs. 

The  overall  scheme  used  to  predict  CZ  range  and  width  is 
to  trace  a series  of  rays  starting  at  the  SLD  with  a zero  de- 
pression angle  ray.  That  first  ray  yields  rQ  which  is  stored. 
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Then  an  iterative  process  is  begun  in  which  the  angle  is  in- 
cremented and  each  succeeding  cycle  distance  determined  is 

compared  to  the  previous  one  until  r . and  9 are  found. 

min  rmin 

(®rmin  stored  use  in  the  cz  9ain  and  TL  program,  to 

be  discussed  later.)  Corrections  are  then  made  to  r . and 

min 

rQ  to  account  for  surface  duct  effects  (if  any)  and  source 
and  receiver  depth  separation  from  the  SLD.  Range  to  the 
inner  edge  of  the  CZ  is  rm^n  plus  corrections , and  range  to 
the  outer  edge  is  Tq  plus  corrections. 

In  the  first  attempt  to  produce  a calculator  program, 
the  cycle  distance  of  ray  #4  (the  ray  just  grazing  the  bot- 
tom) was  compared  to  rQ . The  greater  of  the  two  was  picked 
as  the  basic  distance  for  determining  range  to  the  outer 
edge  of  a CZ.  Later  on,  this  portion  of  the  program  had  to 
be  deleted  to  save  program  steps.  The  final  programs  de- 
veloped ignore  bottom  depth  and  do  not  include  rays  outside 
the  reswept  zone  in  determining  annular  width.  This  is  prob- 
ably a shortcoming  of  the  programs  but  the  seriousness  of 
the  errors  it  causes  will  not  be  known  without  further  study. 

A commonly  applied  rule-of-thumb  states  there  must  be 
a minimum  300  fathoms  of  depth  excess  (water  column  below 
the  DSC)  in  order  to  have  "reliable"  CZ  conditions.  It  was 
observed  that  a fully  developed  reswept  zone  existed  in 
every  case  in  the  locations  studied,  and  separate  calcula- 
tions showed  that  somewhat  less  than  300  fathoms  depth  ex- 
cess was  required  to  complete  the  zone.  Therefore,  a program 
user  should  consider  the  300  fathom  rule-of-thumb  before 
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running  the  range  and  width  program.  With  less  than  300 
fathoms  depth  excess,  the  possibility  exists  for  an  "atypical" 
CS  propagation  situation  where  the  reswept  zone  is  reduced  in 
width  due  to  bottom  ray  limiting. 

Another  program  shortcoming  involves  an  assumption  that 
both  source  and  receiver  would  be  more  shallow  than  the  sec- 
ond DSC  SVP  point  chosen  (depth  D2) . In  other  words,  the 
programs  were  designed  to  allow  for  source/receiver  depths 
within  the  mixed  layer  or  the  first  isogradient  layer  below 
the  SLD.  After  the  five  point  SVP  is  entered  and  the  grad- 
ients computed,  source  and  receiver  depths  are  entered  and 
converted  to  velocities.  The  programs  determine  these  veloc- 
ities (Cg  and  CR)  by  subtracting  an  appropriate  amount  from 
the  velocity  at  the  SLD.  The  amount  subtracted  is  determined 
by  depth  separation  from  the  SLD  and  by  the  gradient  in  either 
the  ML  or  the  first  layer  below  the  SLD.  If  source  or  re- 
ceiver depth  is  greater  than  D^ , sound  velocity  should  be  de- 
termined by  correcting  C2  (the  velocity  at  D2)  and  by  using 
g2  (the  second  layer  gradient).  Since  this  is  not  done,  ve- 
locities for  source/receiver  depths  below  d2  will  be  in  error 
(usually  too  low) . Source  and  receiver  velocity  errors  are 
carried  over  into  Arg  and  ArR  range  corrections.  If  the 
velocities  are  too  low,  the  range  corrections  will  be  too 
large.  This  is  normally  a rather  insignificant  source  of 
total  range  error,  however,  since  Ars  and  ArR  errors  will  be 
a small  fraction  of  the  magnitude  of  those  terms  and  because 
the  range  correction  terms  are  small  to  begin  with. 
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The  mathematics  of  ray  tracing  in  isogradient  layers  is 
quite  straightforward.  By  Snell's  Law, 
c C 

1 m 2 3 LA  constant 

cos  9^  ~ cos  S2  for  each  ray)  1 

the  angle  of  a ray  departing  a layer  can  be  determined  from 
the  angle  of  entry  into  that  layer.  In  Eq.  (2) , is  the 
sound  velocity  where  the  sound  ray  enters  a layer,  8^  is  the 
angle  of  entry,  C2  is  the  sound  velocity  where  the  ray  de- 
parts the  layer,  and  02  is  the  angle  of  departure. 

Rays  travel  in  circular  arcs  within  constant  gradient 
layers,  and  the  radius  of  curvature  is: 


and  9^  are  as  defined  above  and  g^  is  the  gradient  within 
the  layer  (in  this  case,  layer  1)  . 

Finally,  the  horizontal  distance  traveled  by  a ray 
while  traversing  a layer  is: 

Ar  = R (sin  02  - sin  0^)|  (4) 

Figure  5 demonstrates  an  example  application  of  Eqs . 

2 through  4.  It  should  be  noted  that  absolute  value  signs 
are  used  in  Eqs.  3 and  4 because  the  gradient  in  Eq.  3 and 
the  difference  of  sines  in  Eq.  4 may  be  positive  or  negative, 
while  R and  Ar  are  always  positive. 

The  programs  use  these  equations  to  compute  the  horizon- 
tal range  increments  each  ray  accumulates  within  the  four 
layers,  doubles  each  term  Cto  account  for  the  downward  and 
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Figure  5.  Horizontal  Distance  Traveled  Within 

an  Isogradiert  Layer. 
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upward  passes  through  each  layer)  , and  then  sur.is  the  terms 
to  obtain  cycle  distances.  The  fourth  layer  requires  a 
slightly  different  treatment  because  the  rays  become  horizon- 
tal and  then  turn  back  upward  within  that  layer.  Essentially 
the  same  formulas  are  used,  however.  The  equations  are  also 
used  to  compute  the  range  correction  terms. 

In  considering  the  various  possible  ray  paths  between 
source  and  receiver,  it  was  decided  there  were  four  basic 
situations  which  could  occur: 

1)  No  mixed  layer,  both  source  and  receiver 
below  the  SLD.  CDeep/Deep) 

2)  Mixed  layer  present,  both  source  and  receiver 
below  the  SLD.  (Deep/Deep/ML) 

3)  Mixed  layer  present,  both  source  and  receiver 
within  the  layer.  (Shal/Shal) 

4)  Mixed  layer  present,  source  or  receiver  above 
the  SLD,  the  other  below.  (Crosslayer) 

The  only  difference  between  the  first  two  cases  is  the  mixed 

layer  effect  in  case  2.  That  effect  causes  a widening  of 

annuli  due  to  spreading  of  sound  rays  as  they  travel  up  to 

the  surface  and  back  down  to  the  SLD  within  the  layer.  The 

mixed  layer  effect  is  also  included  in  the  third  and  fourth 

cases  above. 

It  was  originally  intended  to  include  all  four  cases  in 
one  range  prediction  program.  Again  due  to  calculator  limi- 
tations, it  was  necessary  to  use  two  programs  to  cover  the 
four  possibilities.  The  first  range  program  (labeled  Deep/ 
Deep)  is  for  cases  1)  and  2)  above  when  both  source  and  re- 
ceiver are  below  the  SLD  whether  or  not  an  ML  exists.  The 
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second  range  program  (labeled  Shal/Shal  or  Crosslayer}  is 
for  use  in  cases  3)  and  4)  above  when  source  or  receiver  or 
both  are  above  the  SLD. 


Formulas  used  to  determine  range  to  inner  edge  of  the 
first  CZ  (RCZi)  and  range  to  the  outer  edge  of  the  first  CZ 
(RCZo)  follow: 

Deep/Deep 

Shal/Shal 

Crosslayer 

Deep/Deep 
Deep/Deep/ML 

Shal/Shal 
Crosslayer 

In  these  equations,  r . and  rA  have  been  previously  defined, 

min  u 

Arg  and  ArR  are  the  respective  horizontal  range  corrections 
which  account  for  source  and  receiver  depth  separation  from 
the  SLD,  and  2 ArQ  is  the  correction  for  mixed  layer  effect. 
Figures  6(a)  through  6(f)  (not  to  scale)  depict  the  RCZi  and 
RCZo  formulas  in  graphic  form.  Ranges  to  second  and  subse- 
quent CZ  annuli  are  taken  to  be  integer  multiples  of  the 
ranges  produced. 

Two  items  of  interest,  both  evident  in  Figs.  6(a)  - 6(f), 
are  worth  mentioning  at  this  point.  First,  acoustical  reci- 
procity is  envoked  and  the  more  shallow  of  source  and  re- 
ceiver is  always  treated  as  "source"  of  the  sound  rays  within 
the  calculator  programs.  Secondly,  only  those  sound  rays 


RCZi 


RCZo 


rmin  ' ArS  ' ArR 


rmin  * 4rS  + “r 


r . + Ar  - Ar„ 

min  s R 


r0  * (2  «„>  + Ars  + irR 


ro  ♦ 2 4ro  + 4rs  ' irR 


r0  + 2 4r0  * ArS  + arR 
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which  experience  no  more  than  one  ocean  surface  reflection 
between  source  and  receiver  are  considered  in  this  model. 

Both  of  these  conventions  are  commonly  applied  to  ray  trac- 
ing models.  Although  they  theoretically  have  little  or  no 
effect  on  model  results,  they  greatly  simplify  the  work  of 
programming  a ray  tracing  model. 

C.  CZ  GAIN  AND  TRANSMISSION  LOSS  MODEL 

In  general,  transmission  loss  is  defined  as  ten  times 
the  logarithm  of  the  ratio  of  sound  intensities  measured  at 
one  meter  from  a source  and  at  range  r from  that  source . 

TL  =*  10  log  ~ 
r 

Intensity  has  units  of  power  per  unit  area.  The  change  in 
intensity  between  one  meter  and  range  r is  due  to  geometric 
spreading  of  the  power  over  a different  amount  of  area  and 
due  to  attenuation  of  some  of  the  power  through  absorption, 
scattering,  diffusion,  etc. 

In  ray  tracing  theory,  it  is  assumed  there  is  no  sound 
power  transfer  across  sound  rays.  Therefore,  the  power  flow- 
ing from  a source  between  two  sound  rays  remains  between 
those  rays  and  travels  out  in  a direction  parallel  to  the 
ray  paths.  Determining  the  portion  of  transmission  loss  due 
to  geometric  spreading  (TL^)  under  this  assumption  reduces  to 
finding  ten  times  the  logarithm  of  the  ratio  of  areas  (at 
range  r and  at  one  meter)  penetrated  by  the  power  between 
the  two  rays  perpendicular  to  the  direction  of  travel. 
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TLg  " 10  io*  A^ 

A mathematical  development  of  this  technique  is  contained  on 
pages  119-121  of  Ref.  1. 

Figure  7 shows  how  this  method  was  adapted  for  use  in 
the  CZ  Gain  and  Transmission  Loss  portion  of  the  calculator 
model  developed.  The  area  (A^)  at  one  meter  from  the  source 
is  the  product  of  area  height  and  area  circumference.  The 
sound  rays  bounding  the  area  above  and  below  are  the  minimum 
and  maximum  departure  angle  rays  which  produce  the  reswept 
zone  in  the  CZ  annulus . The  angular  spread  of  those  rays 
(A8)  in  radian  units  times  the  sphere  radius  (1  meter)  is 
the  area  height.  Cosine  of  the  average  angle  of  departure 
of  the  rays  (8^)  times  the  sphere  radius  times  2tt  is  circum- 
ference of  the  area.  Therefore: 

■ 2ir  A8  cos  9 ^ (m2) 

In  the  CZ,  the  area  (A2)  over  which  the  same  power  is 
distributed  is  also  found  by  a product  of  area  circumference 
and  area  width.  Circumference  is  2*  times  range  to  the  CZ 
(RCZi) . Width  of  the  area  perpendicular  to  the  sound  rays 
is  the  product  of  CZ  annulus  width  (CZW)  and  the  sine  of  the 
average  angle  of  arrival  of  the  rays  at  the  receiver  depth 
(92) . Therefore: 

A2  ■ 2 it  RCZi  CZW  sin  92  (m2) 
and  the  geometric  TL  expression  becomes: 
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RCZi  CZW  sin  9, 

TLg  * 10  lo«  "AT  cos 

Substituting  this  expression  back  into  Eq.  1,  which  is 
the  definition  of  CZ  gain,  and  reducing  to  simplest  form 
yields  the  algorithm  used  to  determine  G in  the  calculator 
model: 

RCZi  A9  cos  0.  ... 

G - 10  lo-J  cZw—  in-T2  ' ’ 

To  implement  this  algorithm,  the  program  has  only  to 
determine  the  angular  terms  since  RCZi  and  CZW  are  available 
from  the  range  program  results.  After  one  of  the  range  pro- 
grams has  been  run,  the  user  loads  the  G and  TL  program  into 
calculator  memory  without  altering  the  contents  of  the  data 
storage  registers  left  from  the  range  program.  Then  the 
iterative  ray  tracing  process  begun  in  the  range  program  is 
continued  in  the  gain  program.  The  angle  of  departure  of 

sound  rays  from  the  SLD  is  incremented  bevond  9 (left  in 

rmin 

storage  from  the  range  program)  and  cycle  distances  produced 

are  checked  for  approximate  equality  with  rQ . In  this  way, 

the  ray  which  completes  the  reswept  rone  is  found,  and  its 

angle  of  departure  from  the  SLD  is  9 . Then  9„  and 

rswp  rswp 

zero  degrees  (the  angle  for  the  ray  producing  cycle  distance 
rQ)  are  converted  to  angles  of  departure  from  the  source 
depth,  0SR  and  9SQ  respectively,  and  angles  of  arrival  at 
the  receiver  depth,  9RR  and  0RQ  respectively,  using  Snell's 
law.  The  angular  terms  in  Eq.  5 are  then  computed  using 
the  following  formulas: 
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9sr 

+ 

9S0 

Deep 

"source" 

9sr 

- 

9 so 

Shal 

"source" 

9sr 

- 

9 SO 

Deep 

"source" 

9sr 

T 

+ 

9 so 

Shal 

"source" 

T 

9ro 

+ 

9rr 

Recall  that  "source"  in  the  model  refers  to  the  more  shallow 
of  source  and  receiver.  Therefore,  the  deep  "source"  forms 
of  these  formulas  are  used  only  after  using  the  Deep/Deep 
range  program.  In  all  other  cases,  the  "source"  is  consid- 
ered to  be  shallow.  Figure  8 depicts  these  angular  relation- 
ships for  the  various  depth  conditions. 

It  should  be  pointed  out  there  are  two  inherent  errors 
in  the  angular  quantities  determined.  First,  the  possible 
source  and  receiver  sound  velocity  errors  mentioned  earlier 
could  cause  the  gain  algorithm  angles  to  be  slightly  off. 

This  would  only  occur  if  depths  greater  than  D.,  were  entered 
for  source  or  receiver  or  both.  Secondly,  6 _ is  found 

for  the  ray  which  has  cycle  distance  equal  to  rQ  at  the  SLD. 
Since  the  actual  CZ  ray  bundle  departs  from  the  source  depth 
(vice  SLD)  and  arrives  at  the  receiver  depth  (vice  SLD) , the 
ray  which  completes  the  reswept  zone  will  probably  be  differ- 
ent than  the  ray  used  and  it  will  have  a slightly  different 
angle  crossing  the  SLD.  These  angular  errors  will  cause  the 
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source 


Shallow  "Source 


Shallow  receiver 


RO  / T2  / \ RR  Deep  receiver 
/ / depth 


Figure  8.  Determining  angular  terms  for 

CZ  Gain  algorithm. 


greatest  CZ  gain  error  in  the  sin  92  term  of  the  algorithm. 
Since  sine  is  directly  proportional  to  angle  at  small  angles, 
an  error  of  a factor  of  two  in  82  (a  quite  possible  event) 
could  cause  a gain  error  of  approximately  3 dB. 

Once  CZ  gain  is  computed  and  stored,  the  sound  frequency 
of  interest  is  entered,  and  the  attenuation  coefficient  is 
calculated  using  Thorpe's  equation  (p.  102,  Ref.  1) : 

(dB/m)  (6) 

In  Eq.  6,  f is  in  kHz,  and  the  constant  in  front  of  the  ex- 
pression converts  attenuation  coefficient  from  dB/kyd  to 
dB/m.  The  program  user  enters  frequency  in  Hz,  and  the  pro- 
gram performs  the  conversion  to  kHz . 

Finally,  the  transmission  loss  in  the  n CZ  annulus 
for  the  frequency  of  interest  (TLfl)  is  determined  by: 

TL  = 20  log  (n  RCZi)  + z (n  RCZi)  - G (7) 

n s 

In  this  equation,  the  subscript  n denotes  the  n CZ  annulus, 
the  range  to  which  is  n times  RCZi. 

After  a range  program  is  run,  and  after  the  gain  portion 
of  the  G and  TLn  program  has  been  completed,  TLn  values  for 
a variety  of  frequencies  and  CZ  annuli  may  be  rapidly  ob- 
tained for  the  SVP,  source  depth,  and  receiver  depth  condi- 
tions entered.  If,  however,  a different  set  of  source/receiver 
depth  conditions  are  also  of  interest,  the  entire  procedure 
beginning  with  the  appropriate  range  program  must  be  per- 
formed again. 


(0.001094) 
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D.  CALCULATOR  PREDICTIONS  COMPARED  TO  I CAPS 

1.  Choosing  the  SVP  Points  for  the  Program 

The  five  point  SVP  limitation  of  the  ray  tracing 
procedure  is  a rather  serious  handicap  in  many  situations. 
Actual  sound  velocity  profiles  not  only  are  curvilinear  in 
overall  shape  but  also  have  many  small  scale  features  and 
they  are  time  varying  functions . Approximating  these  curves 
with  only  four  straight  line  segments  presents  a difficult 
challenge . 

In  general,  matching  the  gradients,  sound  velocities, 
and  associated  depths  are  all  important  in  choosing  SVP 
points.  The  greatest  potential  for  causing  large  range  pre- 
diction errors  occurs  when  the  SVP  contains  an  extensive 
near  surface  layer  with  very  slight  velocity  gradient.  The 
horizontal  distance  traveled  by  a shallow  depression  angle 
ray  within  such  a layer  varies  considerably  with  small 
changes  in  the  gradient  or  layer  thickness.  Under  such  con- 
ditions, then,  it  is  extremely  important  to  match  those 
characteristics  as  closely  as  possible. 

Another  important  item  to  carefully  match  is  the 
sound  velocity  at  the  DSC  axis.  This  velocity  determines 
the  maximum  angle  of  depression  for  each  ray  prior  to  com- 
mencing upward  refraction.  The  horizontal  distance  traveled 
by  a ray  below  the  axis  is  highly  dependent  on  that  angle. 

Table  IV  contains  the  five  point  sound  velocity 
profiles  picked  by  the  author  for  use  in  comparing  the  pro- 
gram performance  to  ICAPS  predictions . Depths  in  the  table 
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5 Points 

FEB 

MAY 

AUG  1 

| NOV 

575 

1524.0 

1225 

1488.5 


MEDITERRANEAN 


10 

1526.6 

30 

1518.0 
100 

1513.0 
800 

1521.5 


; 0 
1508.6 

30 

1497.7 

100 

1483.0 

100 

1484.0 

600 

1475.5  „ 

600 

1476.0 

1750 

1486.5 

1400 

1483.0 

3120 

2460 

0 

1541.1 

0 

1535.5 

125 

1522.5 

125 

1522.5 

550 

1524.5 

625 

1524.5 

1200 

1487.5 

1175 

1489.0 

4760 

4460 

0 

1537.6 

10 

1528.7 

50 

1517.5 

50 

1517.5 

100 

1513.0 

125 

1512.0 

700 

1519.5 

900 

1522.5 

1800 

1290 

Table  IV.  Five  Point  Sound  Velocity  Profiles. 

(Depths  in  meters,  velocities  in  m/sec) 


are  in  meters,  and  sound  velocities  are  in  meters  per  second. 
The  reader  may  want  to  plot  these  points  on  the  graphs  of 
Figs.  1(a)  through  1(1)  so  he  may  see  how  the  four  isogradi- 
ent layers  picked  match  the  ICAPS  profiles.  It  should  be 
pointed  out  that  only  the  initial  selection  of  SVP  points 
was  used  in  the  subsequent  comparisons  of  calculator  model 
results  to  ICAPS  predictions.  Since  an  ASW  aircrewman  using 
the  programs  in  attempting  an  in  situ  prediction  of  acoustic 
conditions  would  not  be  able  to  judge  whether  SVP  point  ad- 
justments would  improve  or  degrade  prediction  accuracy,  it 
was  felt  that  comparing  resuts  of  the  initial  SVP  point  se- 
lection with  ICAPS  would  be  more  meaningful  to  the  objective 
of  developing  the  calculator  programs. 

Comparing  calculator  model  predictions  to  ICAPS  pre- 
dictions in  a definitive  statistical  manner  was  not  done. 

The  main  reason  for  this  was  alluded  to  in  the  preceding  para' 
graphs.  Since  the  SVP  points  entered  in  the  calculator  pro- 
gram must  be  picked  subjectively  by  the  person  using  the 
program  and  since  it  is  unlikely  different  people  would  pick 
the  exact  same  points  off  any  given  SVP,  it  is  clear  that 
calculator  results  can  be  expected  to  vary  from  operator  to 
operator. 

2 . C2  Range  and  Annulus  Width  Comparisons 

The  calculator  range  programs  produce  one  value  each 


for  RCZi  and  RCZo  for  any  given  SVP,  source  depth,  and  re- 
ceiver depth  situation.  Under  the  same  set  of  conditions, 
ICAPS  yields  four  sets  of  RCZi  and  CZW  values,  one  set  for 


each  of  the  four  frequencies  entered.  In  order  to  compare 
the  calculator  performance  to  ICAPS  it  was  first  necessary 
to  reduce  the  ICAPS  predictions  to  one  value  each  for  RCZi 
and  CZW  for  each  SVP/source/receiver  condition.  This  was 
done  by  simple  averaging  to  eliminate  the  frequency  variable 
from  the  ICAPS  range  and  width  predictions . 

Figures  9(a)  - 9(c)  display  range  and  width  compari- 
sons in  graphical  form  for  the  Pacific,  Atlantic,  and  Medi- 
terranean locations  respectively.  In  each  figure  the  double 
barred  lines  represent  the  ICAPS  first  CZ  annuli  predictions 
(averaged  over  frequency) , and  the  single  barred  lines  repre- 
sent the  calculator  predictions.  Numerical  values  for  inner 
and  outer  first  CZ  ranges  may  be  obtained  from  the  scales  at 
the  tops  of  the  figures. 

In  all,  there  were  32  cases  where  these  graphical 
comparisons  could  be  made.  The  following  comments  pertain 
to  those  comparisons : 

a)  In  30  of  the  32  cases  the  calculator  annuli 
overlap  at  least  a portion  of  the  ICAPS  annuli. 

b)  In  14  of  the  32  cases  the  calculator  annuli 
are  completely  contained  within  the  limits  of  the  ICAPS 
annuli. 

c)  In  all  32  cases  RCZi  ranges  predicted  by 
the  calculator  were  greater  than  those  predicted  by  ICAPS. 

In  the  12  Pacific  cases,  the  calculator  RCZi  values  were  ap- 
proximately 1.7  nm  greater  than  ICAPS  on  the  average.  In 
the  12  Atlantic  cases,  the  average  difference  was  approximately 
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Figure  9(a).  Comparisons  of  ICAPS  and  Calculator  CZ 
Annuli  Predictions.  Pacific  Ocean. 
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Figure  9(b).  Comparisons  of  ICAPS  and  Calculator  CZ 
Annuli  Predictions.  Atlantic  Ocean. 
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Figure  9(c).  Comparisons  of  ICAPS  and  Calculator  CZ 
Annuli  Predictions.  Mediterranean  Sea. 
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2.8  nm.  And  in  the  eight  Mediterranean  cases,  1.5  nm  was 
the  mean  difference. 

d)  In  27  of  the  32  cases  the  CZ  width  predic- 
tions from  the  calculator  were  more  narrow  than  the  ICAPS 
predicted  widths.  Three  of  the  five  cases  where  calculator 
CZW  exceeded  ICAPS  CZW  were  from  the  February  SVP  in  the  At- 
lantic location.  That  SVP  contained  a very  deep  (500  meter) , 
nearly  isovelocity  layer  near  the  surface.  In  such  a pro- 
file, CZ  refraction  produces  ray  paths  that  are  spread  over 
a very  wide  (in  this  case  16-TO  nm)  annulus.  Only  the  rays 
which  return  to  the  SLD  within  -he  first  few  nm  at  the  inner 
edge  of  that  annulus  experience  sufficient  convergence  to 
produce  detectable  CZ  gain,  however.  Going  from  inner  to 
outer  edge  of  such  an  annulus  the  CZ  refracted  rays  rapidly 
fan  out  experiencing  progressively  less  convergence  and  pro- 
ducing progressively  less  CZ  gain.  Additionally,  if  the 
bottom  grazing  ray  were  considered,  it  would  be  seen  to  limit 
the  reswept  region  of  this  type  annulus  to  something  far  less 
than  that  indicated.  Since  the  calculator  model  fails  to  ac- 
count for  either  of  these  factors,  it  fails  rather  dramat- 
ically to  produce  a "practical"  CZ  annular  width  from  this 
SVP  type . 

In  summary,  the  calculator  model  produces  CZ  annuli 
that  roughly  agree  with  those  produced  by  ICAPS  in  all  three 
ocean  basins  considered.  Calculator  RCZi  ranges  are  5-10% 
greater  on  the  average  than  the  ICAPS  values.  Calculator 
CZW  predictions  (excluding  the  Atlantic  February  SVP)  are 


40-50%  narrower  than  ICAPS  widths  on  the  average.  And,  the 
Atlantic  February  case  indicates  there  is  at  least  one  SVP 
type  in  which  the  calculator  model  fails  to  produce  even 
marginally  acceptable  results  for  CZW. 

3 . CZ  Gain  Comparisons 

As  with  CZ  range  and  width  comparisons,  it  was  nec- 
essary to  average  the  ICAPS  gain  data  with  respect  to  fre- 
quency before  calculator  gain  predictions  could  be  compared. 
The  estimated  ICAPS  gain  values  in  Tables  II  and  III  were 
thus  reduced  to  one  number  for  each  SVP,  source,  and  receiver 
condition.  Table  V,  CZ  Gain  Prediction  Comparisons,  contains 
numbers  that  represent  the  difference  between  calculator  gain 
predictions  and  the  averaged  ICAPS  values.  Minus  signs  in 
the  table  indicate  those  cases  where  calculator  gain  was  less 
than  the  ICAPS  value. 

As  with  the  range  and  width  comparisons,  the  worst 
agreement  occurred  in  the  Atlantic  winter  SVP  case.  Since 
CZW  is  a term  in  the  gain  algorithm,  the  extremely  wide  an- 
nuli predicted  by  the  calculator  caused  gain  values  to  be 
far  too  low  for  the  three  source/receiver  conditions  assoc- 
iated with  that  SVP. 

Excluding  the  Atlantic  winter  SVP  case,  the  follow- 
ing comments  can  be  made  concerning  the  other  29  CZ  gain 
comparisons : 

a)  Calculator  gain  values  ranged  from  7.3  dB 
lower  to  5 dB  higher  than  the  averaged  ICAPS  values. 
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GICAPS  / ^Gcalc  " GICAPS^ 


Source/ 


SVP 

Profile 

Receiver 

(ft) 

Pacific 

Atlantic 

Mediterranean 

60/60 

13.4/-5.4 

15.5/-10.5 

FEB 

60/400 

11.4/-2.4 

13.0/-7.0 

300/400 

11.2/-2.2 

17.5/-U.5 

60/60 

16.4/  0.6 

I6.5/-O.5 

14.3/  3.7 

MAY 

60/400 

11 .6/-2.6 

14.0/-2.0 

10.7/-1-7 

300/400 

12.7/-0.? 

15.O/-3.O 

I8.3/-6. 3 

60/60 

15-8/  2.2 

16.0/  2.0 

17.7/-0.7 

AUG 

60/400 

12.8/-1 .8 

11.5/  0.5 

10.0/  0.0 

300/400 

12.4/  0.6 

14.5/-1 .5 

60/60 

14. 0/-4.0 

12.0/  5.0 

13.0/  4.0 

NOV 

6c/400 

11.3/-7.3 

11.0/  1.0 

11.0/-3.0 

300/400 

10.7/  1.3 

12.5/  0.5 

13.7/-2.7 

Table  V.  CZ  Gain  Prediction  Comparisons. 


b)  In  22  of  the  29  comparisons,  calculator 
values  were  within  3 dB  of  ICAPS. 

c)  In  nine  of  the  29  comparisons  calculator 
values  were  within  one  dB  of  ICAPS. 

d)  On  the  average,  calculator  gain  values  were 
approximately  one  dB  less  than  ICAPS.  This  result  is  incon- 
sistent with  calculator  CZW  results  in  light  of  the  gain  mod- 
el used.  Since  calculator  CZW  values  averaged  only  slightly 
more  than  half  the  ICAPS  widths,  it  would  have  been  more  con- 
sistent if  calculator  gain  values  turned  out  two  to  three  dB 
higher  than  ICAPS  (.acoustic  power  being  spread  over  less  area 
in  the  CZ  annuli,  other  things  being  equal) . Perhaps  an  ex- 
planation for  this  apparent  discrepancy  is  that  the  calcula- 
tor model  does  not  consider  the  contribution  of  surface 
reflected  energy  adding  to  the  energy  from  upward  traveling 
sound  rays  at  the  receiver  depth.  In  an  actual  CZ  annulus 
the  downward  traveling,  surface  reflected  energy  adds  approx- 
imately three  dB  to  the  CZ  gain  over  much  of  the  annulus 
width.  Apparently,  the  FACT  model  in  the  ICAPS  system  in- 
cludes this  consideration.  It  is  also  apparent  that  neglect- 
ing surface  reflected  energy  in  the  calculator  gain  model  has 
the  effect  of  canceling  errors  that  should  result  from  CZW 
values  being  too  narrow. 

In  summary  of  the  gain  results,  it  can  be  said  that 
the  ray  tracing  technique  used  in  the  calculator  model  worked 
reasonably  well.  Since  three  fourths  of  the  comparison  cases 
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resulted  in  gain  values  within  three  dB  of  the  estimated 
ICAPS  figures,  TL  values  from  the  calculator  displayed  the 
same  close  agreement. 
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IV.  CONCLUSIONS 


A.  LIMITATIONS  OF  THE  MODEL  DEVELOPED 

The  HP-67/97  calculators  used  in  programming  the  CZ  pre- 
diction model  were  stretched  to  their  limits  in  both  data 
storage  and  program  step  capacity.  Although  not  known  for 
certain,  the  author  feels  significantly  more  accurate  results 
would  be  possible  from  a calculator  with  only  moderately 
larger  storage  capacity. 

The  data  storage  limitation  which  allowed  only  five  SVP 
points  to  be  entered  is  quite  restrictive  and  no  doubt  plays 
a large  role  in  the  CZ  range  and  width  inaccuracies  obtained. 

Program  step  capacity  forced  several  short  cuts  to  be 
taken  which  again  would  not  have  been  necessary  with  a mod- 
erately larger  program  memory.  Two  separate  range  and  width 
programs  were  required  due  to  insufficient  program  space  to 
incorporate  tests  for  different  source  and  receiver  depth 
cases.  Also,  source  and  receiver  depths  are  strictly  allowed 
only  within  the  upper  two  SVP  isogradient  layers  because  pro- 
gram space  was  not  available  to  check  for  the  correct  layer 
if  all  depths  were  allowed.  Additionally,  and  perhaps  the 
greatest  source  of  CZW  errors  observed,  program  step  limita- 
tion prevented  incorporating  a method  of  considering  the 
bottom  limited  CZ  sound  ray  in  determining  the  range  to  the 
outer  edge  of  a CZ  annulus.  The  program  developed  ignores 
the  bottom  entirely  and  considers  only  the  reswept  zone  in 
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predicting  annular  width.  Since  calculator  CZW  results  were 
considerably  shorter  than  those  indicated  by  ICAPS,  it  is 
assumed  the  discrepancy  is  due  to  not  considering  CZ  rays 
beyond  the  reswept  zone.  The  first  priority  in  making  im- 
provements to  the  calculator  model,  should  a larger  capacity 
machine  be  implemented,  would  be  incorporating  a better  meth- 
od for  selecting  the  ray  which  defines  the  outer  limit  of 
the  CZ  annulus. 

B.  USEFULNESS  OF  THE  MODEL  DEVELOPED 

The  degree  of  success  in  producing  a useable  CZ  prediction 
model  for  handheld  calculators  must  be  determined  by  consid- 
ering the  objectives  set  forth  in  the  first  section  of  this 
study.  The  central  idea  was  to  ascertain  if  a calculator 
model  would  improve  on  ASRAPS  CZ  prediction  accuracy  in  the 
case  where  BT  conditions  determined  in  situ  differed  from 
those  used  to  generate  -the  ASRAPS  TL  profiles.  Inherent  in 
this  objective  is  the  assumption  that  ASRAPS  TL  profiles 
generated  primarily  from  climatological  data  would  be  in 
error  due  to  lack  of  input  data  accuracy.  Also  inherently 
assumed  is  that  given  identical  input  data  the  calculator 
model  would  produce  less  accurate  results  than  the  digital 
computer  model(FACT)  used  in  A3RAPS(and  ICAPS)  due  to  obvious 
differences  in  data  and  program  capacities.  The  real  question 
then  is  a trade-off  comparison:  Will  the  basically  less  accu- 
rate calculator  model  produce  better  CZ  predictions  with  actual 
environmental  data  than  the  more  sophisticated  digital  com- 
puter model  which  had  only  climatological  input  data? 

L 
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Before  addressing  the  answer  to  this  question,  character- 
istics of  the  calculator  model  developed  will  be  compared  to 
the  list  of  six  desirable  characteristics  described  in  sec- 
tion I. 

1.  Easily  available  input  data. 

The  data  required  are  an  SVP,  assumed  source 

depth,  hydrophone  depth,  and  frequency  of  interest.  The 
only  portion  of  this  information  not  presently  available  to 
ASW  aircrews  is  that  part  of  the  SVP  below  the  1,000  ft 
depth  limit  of  the  AN/SSQ-36  bathythermograph  buoy.  SVP  data 
from  the  surface  to  1,000  ft  (the  area  where  seasonal  and 
diurnal  variations  predominantly  occur)  is  easily  obtained 
from  the  BT  buoy  information. 

2.  Ease  of  program  operation. 

Anyone  familiar  with  HP-67/97  calculator  use 

could  operate  this  program  without  additional  training. 

3.  Output  data. 

The  program  provides  CZ  annulus  range  and  width 
as  well  as  TL  values  for  all  frequencies  of  interest  in  all 
annuli  of  interest. 

4.  Short  run  time. 

To  run  a complete  program  requires  approximately 
10  minutes  once  SVP  data  is  obtained.  Deploying  a BT  buoy 
and  converting  the  temperature  trace  to  an  SVP  would  take  an 
additional  10-15  minutes. 

I 
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5.  Based  entirely  on  acoustic  theory. 


The  program  uses  only  ray  tracing  techniques  in 
producing  its  output  terms. 

6.  Agreement  with  large  computer  models. 

Calculator  CZ  ranges  obtained  averaged  5-10% 
greater  than  ranges  obtained  from  ICAPS.  CZ  width  results 
averaged  only  40-50%  of  those  obtained  from  ICAPS.  And, 
there  was  one  SVP  case  studied  (winter,  Atlantic)  in  which 
the  calculator  CZW  results  were  very  different  from  ICAPS. 
That  SVP  case  was  considered  a failure  of  the  calculator 
model,  and  it  must  be  conceded  the  model  does  not  work  for 
all  CZ  situations.  Excluding  the  obvious  CZW  failure  SVP 
case,  TL  values  from  the  calculator  averaged  about  one  dB 
lower  than  ICAPS  with  extreme  deviations  observed  ranging 
from  -7.3aB  to  +5dB  around  the  ICAPS  values.  Additionally, 
calculator  results  can  be  expected  to  vary  from  operator  to 
operator  since  SVP  points  must  be  picked  subjectively  from 
an  SVP  graph. 

Returning  to  the  main  objective  of  the  study,  the  author 
feels  that  only  half  of  the  trade-off  question  has  been 
answered.  An  easily  operated,  purely  theoretical  model  was 
developed  which  works  for  most  (but  not  all)  CZ  producing 
SVP  conditions.  And  a measure  of  its  accuracy  compared  to 
the  sophisticated  FACT  computer  model  was  obtained.  Yet  to 
be  answered  is  how  inaccurate  ASRAPS  CZ  predictions  are  when 
observed  BT  conditions  differ  from  climatological  conditions. 
This  portion  of  the  question  is  very  difficult  to  answer  and 
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in  fact  would  be  a very  large  study  in  itself.  Generally, 
inaccuracies  must  range  over  a scale  from  insignificant  to 
considerable  as  environmental  deviations  range  from  slight 
to  great.  The  most  likely  variables  affecting  degree  of  in- 
accuracy are  surface  water  temperature,  mixed  layer  depth, 
and  in  layer  and  below  layer  gradients.  The  effects  of  vary- 
ing these  or  other  possible  factors  one  at  a time  or  in 
various  combinations  on  CZ  range,  width,  and  gain  must  be 
known  before  the  entire  question  can  be  answered.  Further, 
definite  magnitudes  of  environmental  factor  deviation  must 
be  determined  so  that  a person  can  judge  when  ASRAPS  inaccu- 
racies are  likely  to  be  greater  than  the  calculator  model 
inaccuracies.  Until  these  points  are  answered  it  would  be 
inappropriate  to  recommend  use  of  the  calculator  model  as  a 
routine  method  of  updating  ASRAPS  CZ  predictions  in  situ. 
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APPENDIX 


HP-67/97  Calculator  Programs  for  Convergence 
Zone  Range , Width  and  Transmission  Loss  Predictions 


Steps  required  to  use  the  programs: 

1.  Deploy  a bathy thermography  buoy  in  the  operating  area 
of  interest. 

2.  Convert  the  BT  buoy  information  to  a sound  velocity 
profile  of  the  upper  1,000  ft  of  the  ocean  area. 

3.  Combine  the  upper  SVP  data  with  a graph  of  climato- 
logical SVP  data  which  depicts  sound  velocity  conditions 
below  the  1,000  ft  level. 

4.  Pick  five  points  from  the  combined  SVP  graph  which 
best  represent  the  deep  sound  channel  portion  of  the  SVP. 

The  first  point  should  be  at  the  sonic  layer  depth,  the  fifth 
point  at  the  bottom  of  the  DSC  where  sound  velocity  equals 
that  at  the  SLD,  and  the  other  three  points  at  points  on  the* 
graph  such  that  when  straight  lines  are  drawn  to  connect  the 
five  points  they  create  a linearly  segmented  SVP  which  matches 
the  actual  SVP  as  closely  as  possible. 

5.  Pick  the  appropriate  CZ  Range  and  Width  program  to  be 
used  as  follows:  If  both  source  and  receiver  are  below  the 
SLD  use  the  Deep/Deep  program.  If  source  or  receiver  or  both 
are  above  the  SLD  use  the  Shal/Shal  or  Crosslayer  program. 

6.  Load  and  run  the  appropriate  Range  and  Width  program 
according  to  the  accompanying  instructions. 

7.  Leaving  the  calculator  power  on  and  data  storage  reg- 
isters unchanged,  load  and  run  the  Gain  and  Transmission  Loss 
program  according  to  its  instructions. 


A word  about  units: 

As  currently  written,  the  programs  use  metric  units; 
meters  for  depths,  and  m/sec  for  sound  velocities.  To  con- 
vert the  programs  for  english  unit  input  data,  feet  for 
depths  and  ft/sec  for  velocities,  the  conversion  factor 
1,852  m/nm  should  be  changed  to  6,075  ft/nm  where  occurring. 
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User  Instructions  for  CZ  Range  Programs 


Step Instructions Input  Keys  Output 


1 

Enter  5-point  SVP  which 

D1 

A 

Dl 

describes  the  Deep  Sound 

X 

Channel.  (D^“SLD , D5» 

ci 

A 

C1 

bottom  of  DSC.) 

D2 

A 

D2 

C2 

A 

C2 

D3 

A 

D3 

C3  • 

A 

C3 

°4 

A 

°4 

C4 

A 

C4 

D5 

A 

D5 

2 

Press  R/S  to  perform  pre- 
liminary calculations 

none 

R/S 

10.0 

3 

Enter  Receiver  Depth  and 

dr 

ENTER 

Source  Depth . 

Ds 

R/S 

RCZi 

4 

Range  to  inner  edge  of 

CZ  (RCZi)  and  range  to 

none 

B 

RCZi 

outer  edge  of  CZ  (RCZo) 
may  be  displayed  by  use 
of  User  Control  keys  B 
and  C respectively. 

Ranges  are  in  nautical 
miles . 

none 

C 

RCZo 
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Storage  Allocation  for  CZ  Range  Programs 
Registers : 


RO; 

D1 

/ 

EAr 

SO; 

gl 

S: 

C.  / cos  9 

0 ' rmm 

Rl; 

C1 

SI; 

g2 

B: 

9SLD  ^ 9 rmin 

R2 

°2 

/ 

CS 

S2: 

g3 

C: 

cos  9 

w 

R3; 

C2 

S3; 

g4 

D; 

COS  9w+l 

R4 : 

°3 

/ 

CR 

S4 : 

E: 

rw 

R5; 

C3 

S5 : 

I: 

Control 

R6  : 

D4 

/ 

rmin 

/ RCZi 

S6: 

R7; 

C4 

S7  ; 

R8  : 

°5 

/ 

ro  / 

RCZo 

S8 : 

R9  : 

g0 

S9  : 

Initial  Flag  Status  and 

Use ; 

0; 

Off, 

Unused 

2: 

Off,  set 

prior  to  rQ 

calculation  in  Deep/Deep 
program  only. 


Is  Off,  On  if  RCVR  is  shallow 
in  Shal/Shal  or  Crosslayer 
program  only. 


3.  Off,  set  by  data  entry 
until  rQ  is  found  in 

S/S  or  Crosslayer  program 
only. 


Display  Status;  DSP  1 


User  Control  Keys; 


A: 

Data  entry 

a; 

B: 

Display 

RCZi 

b. 

C; 

Display 

RCZo 

c : 

D; 

d; 

E: 

e ; 
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Step  Keys  Code 

man 

Stec 

661  $LBLh  21  11 

662  37 Oi  35  45 

663  ISZ1  16  26  46 

694  R-S  51 

Enter  data 

651 

152“ 

•53 

054 

055 

056 

057 
050 

15?“ 

066 

061 

962 

963 

964 
065 
866 
967 
068 
869 

"w!  : =5T 

666  6 66 

66'  2 62 

968  3 TO.’  35  45 

669  ISZI  16  26  46 

Compute 

Gradients 

So 

"373  JTCT3  3T53" 

611  RCLl  36  61 

912  RCL2  36  62 

913  RCL0  36  66 

914  CSB*  23  16  11 

Si 

“373  RIT5  JO  65 
616  RCL3  36  63 
01 T RCL  4 36  94 
916  RCL2  36  62 
019  C SB*  23  16  11 

S2 

“323  3TT7  67 

621  RCL5  36  65 

622  RCL6  36  66 

823  RCL4  36  84 

824  iSBm  23  16  11 

g3 

173“ 

871 

672 

873 

87  4 
875 
976 

177“ 

979 

979 

888 

96 1 

982 

983 

984 

135“ 

986 

987 

988 
969 
996 
691 
992 

SIS  SCE1  3S"3T 

026  RCL7  36  67 

627  RCL8  36  88 

628  PCL6  36  66 

929  CSS*  23  16  11 

«4 

839  CSBk  23  16  12 

Set  1=10 

1 71  SETS  36  66 
932  Xb0*  16-42 
833  CT 06  22  66 

Mixed  layer  ? 

934  ST08  35  68 

935  CTO 1 22  91 

If  no  ML, 

C to  R8 

836  tLBLB  21  M 

•37  RCL9  36  69 

938  * -35 

939  CHS  -22 

949  RCLl  36  81 

94 1 ♦ -55 

•42  STOA  35  11 

c0 

"333  CS6«  2J  ib  13 

944  RCL*  36  11 

945  1 91 

94 6 CSB4  23  16  14 

94 7 $706  35  06 

•40  S7>6  35-55  #6 

ML  present, 

2 ^rQ  to  SS 

993 

994 

895 

896 
697 

998 

999 
100 

Set  1=10 

Range  Program  (Dee' 


86 


36  66 

~45  Compute  C0 
36  45  * a 

-35 
36*  01 
-55 
35  64 


Initialize 
^r  routine 


QSBc  23  16  13 

RCLJ  36  63 

RCLC  36  13  2 ^r, 

CS8d  23  16  1 4 1 

ST *6  35-55  66  I 


Cl 

EHTt  -21 

EHTt  -2i 

RCL  5 36'  05 

RCLO  36  14 
CSB4  23  16  14 
ST+6  35-55  66 
ST *9  35-55  66 


2 ^r. 


2 Ar, 


Step  Key a Code  I Explanation 


191 

ST*8 

35-55  00 

192 

S T*9 

35-55  00 

■ ir—w 

105 

ST03 

72  03 

196 

6T04 

22  04 

~Tsr~ 

n mw 

■ ri  ii 

106 

RCL9 

36  00 

109 

ST+8 

35-55  08 

110 

1 

01 

111 

♦ 

-55 

112 

stoe 

35  15 

wnmniKm 

JTW 

114 

RCL8 

36  00 

115 

RCLE 

36  15 

116 

X6Y? 

16-35 

117 

CTOS 

22  05 

118 

Rl 

-31 

119 

sm 

35  15 

■ F Mil  ■ 

36  *2 

121 

• 

-62 

122 

5 

05 

123 

♦ 

-55 

124 

STO0 

35  12 

125 

COS 

42 

126 

STOC 

35  13 

127 

6702 

22  92 

21  iS~ 

1 29 

ST06 

35  96 

JSTT 

131 

0 

-62 

132 

5 

05 

133 

- 

-45 

134 

COS 

42 

135 

STOA 

35  11 

■ 1 ■ f 1 ■ 

MftM 

137 

RCL2 

36  02 

136 

l 

01 

139 

6S8d 

23  16  14 

149 

ST *9 

35-55  00 

141 

OSZI 

16  25  46 

MUM 

Wild 

23  16  13 

143 

RCL2 

36  02 

144 

PC  Lh 

36  11 

145 

6SB4 

23  16  14 

146 

ST-6 

35-45  06 

147 

0SZ1 

16  25  46 

Set  1=10 


found  ? 

Store  Tq 
<&  first  r 


Stei 

■>  Keys  Code 

151 

Co  6* 

23  16  14 

152 

S'*8 

35-55  08 

153 

OSZI 

16  25  46 

154 

CS6c 

23  16  13 

155 

RCL4 

36  84 

156 

RCLA 

36  11 

157 

C$84 

23  16  14 

158 

ST-6 

35-45  06 

r . found  ? 
min 


otore 
next  rw 


Increment 


36  94 
91 


SSL D * 

COS  &r.Tr 


Next  r routin 


Store  r . 


min 


Initiate 

*%  4 ^a 

corrections 


corrections 
to  r„ 


to  r_ 


^rR 

corrections 


13  1*  12 


36  96 
51 


*LBL& 

21  16  11 

- 

-45 

Ri 

-31 

- 

-45 

Rt 

16-31 

f 

-24 

mi 

35  45 

1SZI 

16  26  46 

tLBLc  21  16  13 
PCL1  36  91 

ENTf  -21 

EHTt  -21 

RTH  24 


Explanation 


to  rr 


to  r . 
min 


Convert 
RCZi  & RCZo 
to  nm 


et  1=10 


Display  RCZi 


Display  RCZo 


Gradients 

Subroutine 


Set  1=10 
Subroutine 


AT 

Initiation 

Subroutine 


Subroutine 


CZ  Range  Fro 


Deep/Deep  Case 


Explanation 


Explana  Lion 


Step  Keys  Code 


Step  Keys  Code 


001 

*lELh 

21  11 

AD*? 

rv4 

$70  i 

35  45 

005 

IS  21 

IS  26  46 

004 

R/S 

51 

w 

• 

000 

0 

00 

00“ 

£ 

02 

000 

STOi 

35  45 

009 

fSJJ 

16  26  46 

010 

RCL3 

36  83 

011 

RCLl 

36  81 

012 

RCL2 

36  82 

015 

RCL8 

36  00 

014 

;spa 

23  16  11 

1W 

931 


139" 


04  J 
04  Z 
943 

J&L 


W 

040 


Enter  data 


Compute 

Gradients 


Sr 


Si 


175 

RCL5 

36 

05 

016 

*ci3 

36 

03 

817 

RCL4 

36 

04 

818 

RCL2 

36 

02 

019 

GSBo 

23  16 

11 

g, 


020 

021 

PCL5 

56  05 

022 

RCL6 

36  06 

023 

RCL4 

36  04 

024 

£S8« 

23  16  11 

S' 


r-w»« 

026 

0CL“ 

36  87 

02" 

ffCLfi 

36  88 

026 

RCL6 

36  86 

029 

CS8* 

23  16  11 

g4 


T 
sroi 


00 
J5  46 


Set  1=9 


nm 

36  00 

933 

RCL9 

36  09 

034 

X 

-35 

035 

CHS 

-22 

036 

RCLl 

36  01 

037 

♦ 

-55 

038 

STC0 

35  11 

Compute  C. 


13 ST 

RCLh 

1 

iSBJ 

ST08 

sm 


045 

W 


SSBb 


mm 

36  11 

01 

23  16  14 
35  08 
JSzSLM 


2 ^rQ  to  R8 


■w 

x:y 


23  10  12 

II 


73=31 

-41 


Set  1=1 0 
Enter  DR  ST 


Reciprocity 
test 


049 

050 


/?CL0 


36  00 
-45 


Compute  C., 

d 


051 

852 

057 

054 

055 


RCL9 

x 

RCLl 

* 

$762 


36  09 
-35 
36  01 
-55 
35  02 


“SFT 

JJT 

05“ 

RCL8 

56  00 

056 

X>Y? 

16-34 

059 

6700 

22  00 

Set  Flag  0 

060 

67Q1 

22  01 

for  sha3.1ow 

061 

*LBL8 

21  00 

Receiver 

062 

$ FI 

16  21  01 

063 

DSZI 

16  25  46 

and 

064 

*LBLl 

21  01 

065 

- 

-45 

Compute  C0 

066 

RCLl 

56  45 

* 

0b » 

X 

-35 

066 

RCLl 

36  01 

069 

♦ 

-55 

$704 

HJ.±j 

871  £SBk>  £3  16  l£ 


Set  1=10 


BhUPH 

B*!  3d 

K 3(8 

■;  nT<1 

H /Cl 

Be  ci 

2 


TTT" 

RCLS 

36  05 

080 

twn 

-21 

061 

£Hff 

-21 

082 

RCL5 

56  05 

085 

36  14 

084 

iySS'd 

23  16  14 

085 

$7*8 

35-55  60 

886 

$7*9 

35-55  08 

15T" 

RCL5 

36  05 

688 

EH7t 

-21 

069 

tNTt 

-21 

090 

RCL7 

36  07 

091 

RCLL 

36  14 

092 

6SB4 

23  16  14 

093 

$7*8 

55-55  80 

094 

$ 7*8 

55-55  00 

f:HJi 

RCLl 

j6  01 

896 

RCLB 

36  14 

887 

* 

-24 

098 

RCL  i 

56  45 

099 

* 

-2* 

180 

RClD 

36  14 

2 ^r-. 


2 &r. 


Ar4 


Step  Keys  Code 

161 

CSfi* 

23  16  15 

162 

X 

-35 

163 

ST+8 

35-55  66 

165  fSBfc  23  16*  12  I Set  1=10 


J 


3 rQ  found  ? 


lie  PCLti  36  66 

111  $l>$  35-55  98 

112  1 81 

113  * -55 


1J4  S TOE 


35  15 


126 

fU 

-31 

121 

STOE 

35  15 

123 

• 

-62 

124 

5 

65 

125 

♦ 

-55 

126 

sroe 

35  12 

12? 

cos 

42 

128 

sroc 

35  13 

12$ 

T7 5F~ 

il“52P 

w*m 

21  65 

131 

ST06 

35  66 

23  16  13 


Store  rQ 
& first  r. 


r . found  ? 
min 


Store 
next  rw 


Increment 
®SLD  4 
005  ®SLD 


ext  r routi 


Store  r . 


Initiate 

^r3  4 &R 
corrections 


corrections 


to  r. 


to  r 

min 


Step  Keys  Code  1 Explanation 


16*  25  46 
23  16  13 
36  64 
61 

23  16  14 
Id  25  46 
16  23  81 
-22 
35-55  66 


23  lt> 
16  25 
16  23 


13 
36  64 
36  11 
23  16  14 
16  23  61 
-22 
35-45  66 


<SrR 

corrections 
to  r„ 


to  r_ 


Convert 
RCZi  & RCZo 
to  nm 


Set  1=10 


isplay  RCZi 


isplay  RCZo 


Gradients 

Subroutine 


Set  1=10 
Subroutine 


Initiation 

Subroutine 


16  23  61 


CZ  Range  Program  (Shal/Shal  & Crosslaver  Ca 


' 


— — — f 


Step  Keva  Code  Explanation 


Keys  -ode  I iixnlaria t j on 


201 

X 

-35 

292 

x:y 

-41 

293 

T 

-24 

294 

STOD 

35  14 

295 

6 SB* 

23  16  15 

296 

RCLC 

36  13 

29? 

es  Be 

23  16  15 

296 

- 

-45 

299 

X 

-35 

216 

RCLC 

36  13 

211 

r 

-24 

212 

RCLi 

36  45 

213 

IS  Z1 

16  26  46 

214 

* 

-24 

215 

hBS 

16  31 

216 

RTN 

24 

21? 

tLBLa 

irTTTT 

216 

X* 

53 

219 

CHS 

-22 

229 

1 

61 

221 

* 

-55 

222 

rx 

54 

223 

RTH 

24 

.ZL<. 

*'S 

IL 

Ar 

Subroutine 


cos  to  sin 
Subroutine 


CZ  Range  Fro 


Shal/Shal  & Crosslayer  Cases' 


User  Instructions  for  CZ  Gain  and  Transmission  Loss  Program 


Step Instructions Input  Keys  Output 

1 Run  one  of  the  CZ  Range 
Programs.  Leave  the 
calculator  on  and  all  of 
the  data  storage  regis- 
ters unchaged. 

2 Load  the  CZ  Gain  and 
Transmission  Loss  Program. 

3 Enter  a zero  if  the  Shal/  0 A G (dB) 

Shal  or  Crosslayer  program 

was  used,  or  enter  a one  or 

if  the  Deep/Deep  program  1 

was  used. 

4 Enter  a frequency  of  in-  f (Hz)  B a (dB/m) 

terest  in  Hz. 

5 Enter  the  number  of  a CZ  n C TLn  (dB) 

annulus  of  interest.  (Step 

5 may  be  repeated  for  as 
many  CZ  annuli  as  desired.) 

6 To  determine  the  TL„ 

n 

values  for  different 
frequencies,  return  to 
step  4. 
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Storage  Allocation  for  CZ  Gain  and  Transmission  Loss  Program 
Registers : 


increments  used  in 
finding  esLD  . 

Display  Status:  DSP  0 

User  Control  Keys: 

A:  Compute  G 
B : Compute  a 

C : Compute  TLr 
D: 

E: 


until  Tq  found. 


n 


a: 
b : 
c : 
d: 
e : 


k 

' 


« 

RO: 

Ar 

SO: 

*1 

A: 

9SO  / 0RO  / f 

Rl: 

C1 

SI: 

CM 

B: 

0SLD  / 9rswp 

• 

R2: 

cs 

S2: 

93 

C: 

cos  0w  / 

R3: 

C2 

S3: 

^4 

cos  9 

rswp 

R4 : 

CR 

S4 : 

a terms 

D: 

cos  9w+l 

R5  : 

C3 

S5  : 

E: 

r0  ^ °SR  / 0RR 

R6  : 

RCZi 

S6  : 

I: 

Control 

R7 : 

C4  / a 

S7 : 

R8: 

RCZo 

S8  : 

I 

R9  : 

G terms  / G 

S9  : 

Initial  Flag  Status  and 

Use : 

• 

0: 

Off,  On  for  shallow 
sound  source. 

2: 

Off, 

unused. 

Is 

Off,  On  to  decrease 

the 

3: 

Off, 

set  by  data  entry 

93 


Keys  Code 


Step  Keys  Code  I gxplana  t.jon 


115  COS 


115 

S IN 

41 

129 

ST*9 

35-24  99 

aim 

RCL8 

36  98 

122 

RCL6 

36  06 

1 23 

STx? 

35-35  99 

STr9  35-24 


126  RCL9 


i3s  *LBLB 


Convert  f 
to  kHz 

Store  f2 

Attenuation 

Coefficient 


.First  term 


Second  term 


Store  a 

16  56  iTcTear  S4) 


Convert  a 
from  dB/kyd 
to  dB/m 


j5 -35  97 


173 

tLBLC 

21  13 

174 

RCL6 

36  96 

175 

A 

-35 

176 

EHT  t 

-21 

Display  a 


nRCZi 

201og(nRCZi) 
a (nRCZi) 


Compute  & 
Display  TL 


tLBLa  21  lb  11 
STOC  35  13 
\ -35 

■\CY  -41 

t -24 

3700  55  14 

QSSb  23  16  12 
RCLC  36  13 
CSBb  23  16  12 
-45 

x -35 

PCLC  36  13 
t -24 


Ar 

Subroutine 


CZ  Gain  and  Transmission  Los 


95 


6s  ru 


r 
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